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Abstract 

We apply a multiple-time version of the reductive perturbation method to 
study long waves as governed by the shallow water wave model equation. As 
a consequence of the requirement of a secularity-free perturbation theory, 
we show that the well known N-soliton dynamics of the shallow water wave 
equation, in the particular case of a = 2f3, can be reduced to the N-soliton 
solution that satisfies simultaneously all equations of the Korteweg-de Vries 
hierarchy. 
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I. INTRODUCTION 



The shallow water wave (SWW) equation, 

u xxxt + au x u xt + (3u t u xx u x t u xx = , (1) 

is known to present solitonic solutions for two families of values of the parameters a and (3. 
The first family, given by 

a = 2(3, (2) 

was shown to be integrable by the inverse scattering method by Ablowitz, Kaup, Newell and 
Segur. 1 ) The second family, given by 

a = P, (3) 

was shown to possess N soliton solutions by Hirota and Satsuma. 2 ^ In addition, it has been 
shown 3 ' 4 ) that Eq.(0) passes the Painleve test if and only if either condition (|2]) or (|3]) is 
satisfied, which is an indication that this equation might be integrable only for these two 
cases. We would like also to remark that Eq.(|l|) may be written as an integrodifferential 
equation: 

U xxt + aUU t -(3U x U t dx'-U t -U x = 0. (4) 
J x 

In this form, the SWW equation was shown to be derivable in the Boussinesq approximation 
from the classical shallow water theory, 4 ) but with the condition 

(3 = 2a, (5) 

which is neither of the above discussed cases. 

In this paper, we are going to apply a multiple-time 5,6 -' reductive perturbation method 7 ^ 
to study the long-wave limit of the SWW equation. This method, when applied to non- 
linear weak dispersive general systems, introduces a connection between the equations of 
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the Korteweg-de Vries (KdV) hierarchy and the elimination of the soliton-related secular- 
ises appearing in the evolution equations of the higher-order terms of the wave-field. As 
a consequence of this approach, we will show that the well known solitary-wave solution to 
Eq. ([!]), existing for any a and (3, can be written as a solitary-wave satisfying simultaneously 
all equations of the KdV hierarchy. This result is a consequence of an intrinsic property of 
the KdV solitary- wave, which allows for a truncation of the perturbative series. Then, we 
will show that, for the case a = 2(3, the N soliton dynamics of the SWW equation can be 
reduced to the dynamics of the N soliton solution of the KdV hierarchy. This result 8 -* comes 
from the fact that the SWW equation with a = 2(3, and the KdV equation are both related 
to the same eigenvalue problem in the AKNS scheme. ^ On the other hand, for a = (3, the 
dynamics underlying the N soliton solutions of the SWW equation has a different nature, 
being similar to the dynamics governing the N soliton solution of the Boussinesq equation. 3 ) 
This case will not be considered in this paper, which is organized as follows. In Section II 
we introduce the multiple-time perturbation scheme based on the long-wave expansion of 
the dispersion relation. In Section III we obtain the first few evolution equations, and in 
Section IV we show how the KdV hierarchy equations appear as a compatibility condition. 
The case of a solitary-wave solution for the KdV hierarchy is considered in Section V, and 
in Section VI we describe how the N soliton solution of the SWW equation with a = 2(3 can 
be obtained as the N soliton solution of the KdV hierarchy. And finally, in Section VII we 
summarize and discuss the results obtained. 



II. LONG WAVE EXPANSION AND THE MULTIPLE TIMES 



The linear dispersion relation of the SWW equation is given by 



To study its long- wave limit, we put 

k = en , (7) 
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with e a small parameter, and we expand the dispersion relation according to 

uj(k) = en - e 3 K 3 + e 5 n 5 - e 7 k 7 H . (8) 

Introducing this expansion into the plane wave solution of the linear SWW equation, we get 

u = a exp % ke{x - t) + n 3 e 3 t - K 5 e 5 t + K 7 e 7 t , (9) 

where a is a constant. Inspired by this solution, we introduce a slow space 

£ = e(x-t), (10) 
as well as an infinite sequence of properly normalized slow time variables 

r 3 = e 3 t ; r 5 = -e 5 t ; r 7 = e 7 t ; etc. (11) 
Accordingly, we have that 

i = 4t (12) 



and 



dt~ e di +e dr 3 dr 5 +€ dr 7 ' {U) 



As we are going to see later, the slow time normalizations introduced in the definitions (]TT[) 
are crucial in the sense that they will allow for an automatic elimination of the secularities 
appearing in the evolution equations of the higher order terms of the wave-field. 9 -* 

III. MULTIPLE TIME EVOLUTION EQUATIONS 

We expand now the wave-field u according to 

u = eu = e(u + eui + t 2 u 2 H J , (14) 

and substitute it, together with Eqs.fll^) and flI3|), in the SWW equation (TJ). The result is: 



9w / <9 2 3 d d 5 d d \ A 

+ "a?(- e ae + ^-^8S + -)« = - (15) 

We proceed then to an order-by-order inspection of this equation. At order e° we get 

dV) . d 2 u , du d 2 u 

-W + d^ + {a + (3) ^W = °- (16) 

We make now, in all components u n of the wave-field expansion, the following transforma- 
tion: 

du n 6 



di a + {3 n ' 
Consequently, Eq.([U|) acquires the form 



n = 0,1,2,... . (17) 



dv d 3 v dv . , 

d7 3 + W~ 6vo ^ = ' (18) 

which is the KdV equation in the time r 3 . 
At order e 1 , Eq.([l5|) gives: 

<9 V . d 2 ui duod 2 U! duid 2 u 

-W + d^ + {a + (3) ^W + {a + ^W = °- (19) 

Transforming Uq and U\ according to Eq . (|TTD , we obtain 



that is, t>i satisfies a homogeneous linearized KdV equation in the time r 3 . 
At order e 2 , Eq.flT5|) gives: 

d A u 2 d 4 UQ d 2 U2 d 2 u 



<9£ 4 <9£ 3 <9r 3 d£dr 3 <9£<9r 5 
u d 2 u 2 du\d 2 u\ du 2 d 2 u \ 

dj~de + ~W~dtF + ~W~d¥ ) ' ~ d£ d£dr 3 ' " <9£ 2 dr 3 



fdu d 2 u 2 duid 2 ui du 2 d 2 u \ du d 2 u d 2 u du _ /01 , 



Transforming uq, U\ and u 2 according to Eq.([T7|), and using the KdV equation flTSP to express 
v 0t3 , we obtain 
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dv 2 d 3 v 2 r d dv 



18a + 24/3' 



+ 



12a + 6(3' 
a + 13 



v 



d 3 v 

w 



a + f3 
36a + 18/3 



dv d 2 v 
) d 5 i 



2 dv d 5 v , 
v o - ^rr + 6wi— . 



(22) 



a + f3 

This equation involves the evolution of vq in the time T5, which is not known up to this 
point. Moreover, when vq is assumed to be a solitary- wave solution of the KdV equation, 
the source term d 5 v /d^ 5 will be a secular producing term. Therefore, before continuing we 
have to solve these two problems. 



IV. THE KDV HIERARCHY AS A COMPATIBILITY CONDITION 

As we have already seen, vq satisfies the KdV equation in the time r 3 . The evolution 
of Vq in the higher order times r 2n+ i can then be obtained in the following way. 5 ) First, to 
have a well ordered perturbative scheme, we impose that each one of these equations be e- 
independent when passing from the slow (t> , £, r 2n+ i) to the laboratory coordinates (v, x, t). 
This will select all possible terms to appear in the equation for v 0T2n+l . Then, by imposing 
the natural (in the multiple time formalism) compatibility condition 

(^0r 3 ) = , (23) 

V ' ^2n+l V ' T3 

it is possible to determine the above constants in terms of a 2n +i, which is left as a free- 
parameter. As it can be easily verified through an explicit calculation, 5 ) the resulting equa- 
tions are exactly those given by the KdV hierarchy. As an example, let us consider the 
evolution of v in time r 5 . From the e-independence requirement, v 0r5 is restricted to be of 
the form 

V0t 5 = "5^0(50 + A>^0& + 75^0^0^ + S 5 VqVq^ , (24) 

where 015, fa, 75 and £5 are constants. The compatibility condition 

(v 0T[i ) = (v 0rr> ) , 
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then, determines /?s, 75 and 5$ in terms of 0:5, which is left as a free-parameter, yielding the 
equation 

vot 5 = fo(5£) - 10t» fo^ - 20t;o^Uo« + 30f \^ • (25) 
Analogously, for the evolution of Vo in T7, we obtain 

vo T7 = - v 0( rt) + 14w fo(55) + 42^f (45) + 140(w ) 3 woc 

+ 70v ^v ^ - 280v v ^v ^ - 700og) 3 - 70v 2 v ^ ■ (26) 

Equations (|25| ) and fl2"6|) are the first two higher-order equations of the KdV hierarchy. The 
same procedure can be used to generate any higher-order equation of the KdV hierarchy. 
The right-hand side of these equations would in principle appear multiplied by the cor- 
responding free-parameter left at each order, which would account for different possible 
slow time normalizations. However, since we have already defined appropriate slow time 
normalizations, these parameters were taken to be 1 in order to have an agreement with 
the normalizations introduced in Eq. (|TT|) . This is an important point since it will allow 
for an automatic elimination 9 ) of the soliton related secular producing terms present in the 
higher-order perturbation equations, which are always of the form 10 -* 

^o(2n+i)£ ; n = 0, 1, 2, . . . . (27) 



V. THE SOLITARY WAVE SOLUTION 

In this section, we are going to consider some specific solutions to our equations. First, 
we will assume for v\ the trivial solution 

«i = 0. (28) 

Then, by taking into account that the evolution of vq in the time T5 is now known to proceed 
according to the fifth-order equation of the KdV hierarchy, given by Eq.(|25|), we obtain for 
Eq.(H 
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2a — 4/5 r 2 1 

V2r 3 + V 2m - 6{V V 2 )Z = —— ~3V V 0( : + V V ^ - Vo^Vott . (29) 

a + p 1 J 

We see in this way that the use of the properly normalized KdV hierarchy equation to 
express vq T5 has automatically canceled out the secular producing term fo(50 from Eq.(^). 
Therefore, besides satisfying the KdV equation in the time t$, the wave-field vq must satisfy 
also the fifth-order KdV hierarchy equation (p5[) in the time r 5 . Actually, in order to 
have a secularity-free perturbation theory up to any higher order, Vq will be required to 
satisfy simultaneously all equations of the KdV hierarchy, each one in a different slow time 
coordinate. 

Next, we assume for vq a KdV solitary-wave solution. As vq must satisfy simultaneously 
all equations of the KdV hierarchy, we take 

v = -2/t 2 sech 2 [k£ - 4k 3 t 3 + 16k 5 t 5 - 64/t 7 r 7 H ] . (30) 

Now, using this solitary-wave solution, one can easily check that the right-hand side of 
Eq.(p9|) vanishes, implying that, along with v±, v 2 also satisfies a homogeneous linearized 
KdV equation: 

V2t 3 + v 2 tft - 6(v v 2 )^ = . (31) 
Like in the case of v\, we will assume for it the trivial solution 

v 2 = . (32) 
We proceed now to the next order, in which case Eq. (|T5D gives 

U3(40 - M3£r 3 + (a + (3) (M05« 35 ) f = , (33) 

where, due to the fact that v\ = v 2 — 0, we have already used that u\ = u 2 = 0. Then, by 
transforming Uq and M3 according to Eq.(^), we get the following equation for V3: 

Vzt Z + V3(tf ~ 6 (vov 3 ) e = . (34) 
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We see in this way that also V3 satisfies a homogeneous linearized KdV equation, and as in 
the previous cases we choose for it the trivial solution 

v 3 = . (35) 

This is actually what will happen at every higher order: after using the equations of the 
KdV hierarchy to express vo T2n+1 , and after assuming the solitary-wave solution ( |3~0|) for t> , 
the evolution of v n in the time T3 will be given by a homogeneous linearized KdV equation. 
Consequently it will always be possible to choose the trivial solution 

v n = ; n= 1,2,3,... . (36) 

It should be mentioned that the evolution of v n (n = 1, 2, 3, . . .) in each higher-order time 
r 2n+3 is always given by a linear partial differential equation, which becomes homogeneous 
after the choice t> n _i = 0. Therefore, the solution ( |3~6"D is in indeed always possible and 
compatible with the whole perturbative scheme. 

We return now to the laboratory coordinates (x, t, u). First, we recall that we have 
assumed for u the expansion given by Eq.(|H|). Moreover, we have obtained a particular 
solution in which u\ = u<i = u 3 = ■ ■ ■ = 0. Therefore, the laboratory wave-field u is given 
simply by 

u = eu . (37) 

On the other hand, Vq was taken to be given by Eq. (p0[) , which is a solitary-wave satisfying 
all equations of the KdV hierarchy. By using the transformation (|17D, we can obtain the 
corresponding uq\ 

12k 2 r™ 



Uq 



a + (3 J$ 

Integrating and using Eq.(^), we get 
12eK 



poo r 1 

J sech 2 - 4/? 3 r 3 + 16k 5 t 5 - 64k 7 t 7 + • - -J d? . (38) 



U — n 

a + (3 



1 - tanh« - 4/t 3 r 3 + 16k 5 t 5 - 64/t 7 r 7 + • • •) . (39) 



We drive now u back to the laboratory by eliminating the slow variables (k, £, T2 n +i), which 
are related to the laboratory ones (k,x,t) according to equations (0), ([Up and ([T1D . The 
result is 

12k 



u 



1 - tanh (kx - k{l + 4k 2 + 16F + Q4k b + • • -}t 



a + (3 

But, the series appearing inside the curly brackets can be summed, yielding finally 



(40) 



12k 



u 



a + p 



1 — tanh kx 



kt 



l-Ak 2 



(41) 



which is the well known solitary-wave solution of the SWW equation ([!]), valid for any value 
of a and (3, except of course a = —(3 . 



VI. TWO OR MORE SOLITON SOLUTIONS 

The elimination of secular producing terms through the multiple-scale method is not 
always possible in the case where a N-soliton solution is assumed for the KdV equation 
because of the possible appearance of obstacles in the perturbative series. 11 ) Obstacles are 
higher-order effects coming from the nonintegrability of the original equation which, when 
present, preclude the existence of uniform asymptotic expansions. An example of such a 
manifestation can be found in Ref. [12]. Here, however, we consider a particular case where 
no obstacles are present, namely the SWW equation with a = 2(3. 

Let us then return to the perturbation theory, but considering now , the case in which 
Vo is assumed to be the two-soliton solution of the KdV equation (|18|), which is given by 

K \ ~ K \ + K \ cos1i(2k 2 72) + k-2 cosh(2/ti7i) 
[(«2 - «i) cos1i(aci7i + K2I2) + («2 + K i) cosh(ACi7! - K2I2)] 2 ] 

with 

Ti = £-4/^T3 ,i = 1,2. (43) 

In this case, the whole perturbative scheme follows closely the one when Vq is a one-soliton 
solution to the KdV equation. In particular, provided the evolution of vq in the time T5 be 

10 



Vq = -4(«| - k\) 



given by the fifth-order KdV hierarchy equation ( p5|) , the equation for v 2 in the time T3 will 
again be given by 

2a — 4/3 r 2 1 

V2t 3 + v 2m - 6{vov 2 )z = —— -3v voz + v vo$tf - v 0( :Vott . (44) 

a + p 1 J 

However, when the two-soliton solution (H2T) is substituted in its right-hand side, it does 



not vanish as it did when vq was assumed to be the solitary-wave solution fl30l) . But, as can 
be easily seen, for the specific case of a = 2/3, Eq.fl4"4j) will be homogeneous independently of 
the solution chosen for vq. Therefore, also for the case in which vq represents a two-soliton 
solution, it is possible to choose the trivial solution v 2 = 0. This is actually what happens 
at any higher order. In fact, to obtain a secularity-free perturbation theory, we first notice 
that must satisfy simultaneously all equations of the KdV hierarchy, which means that 
instead of (|43|) we have (i = 1,2): 

Ti = £ - 4/c?t 3 + 16«?r 5 - 64^r 7 + • • ■ . (45) 

Then, for the case a = 2/3, the evolution of v n in the time r 3 will always be given by a 
homogeneous linearized KdV equation, whose solution we assume to be 

v n = ; n = 1,2,3, ... . 

In other words, the perturbative series truncates, and an exact solution u = e 2 «o can be 
found also for a two-soliton solution. 

The return to the laboratory coordinates (x, t, u) proceeds in the very same way as for 
the case of the solitary-wave. Writing (i = 1, 2) 



Hi — e k.; 



1 1 



and using Eqs.(^) and (|TTD, we find that, in the laboratory coordinates, the two soliton 
solution of the KdV hierarchy is given by 



12 

u = eu = —s(h + k 2 ) 

a + p 



1 (h ~ k 2 ) tanh(fc 2 72) 

k 2 — k\ tanh(A;i7i) tanh(A;272) 
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(46) 



where now 



* =i -ot i = 1 > 2 - < 47 > 



This is exactly the two-soliton solution of the SWW equation when a = 2(3. Therefore, 
like in the case of a solitary-wave, the two-soliton solution of the SWW equation is nothing 
but a two-soliton satisfying simultaneously, in the slow variables, all equations of the KdV 
hierarchy Actually, as a similar procedure shows, this statement holds for the N-soliton 
solution as well. 

The above results indicate that there exists a close relation between the KdV hierarchy 
and the SWW equation with a = 2(3. The root of this relation becomes more explicit in the 
AKNS scheme. x ' Indeed, as a by-product of this scheme, it has been shown that the family 
of equations 

q t + C(4L)q x = 0, (48) 

where 

Id 2 1 r°° 

L = -7^-q + ^ d v, ( 49 ) 



Adx 2 * 2 jx 

and C(k 2 ) = u(k)/k, with u(k) a given dispersion relation, is integrable, having all of them 
the same associated eigenvalue problem. In the specific case of C{k 2 ) = 1/(1 + k 2 ), Eq.(f4"8"D 
is just the SWW equation with a = 2(3. On the other hand, if we expand C{k 2 ) according 
to 

C(k 2 ) = 1 - k 2 + k A - k 6 + ■ ■ ■ , (50) 

we get 

q t + q x - AL(q x ) + lQL 2 (q x ) - ■ ■ ■ = . (51) 

When passing to the multiple-time slow variables (go, £, T 2n+ i), with q = e~ 2 q, one can 
immediately see the emergence of the KdV hierarchy as 

%r 2n+1 = (4L) n g 05 
12 



does generate the hierarchy. 13 ) It should be mentioned that the above result is independent 
of any specific solution to the KdV hierarchy. 

VII. FINAL REMARKS 

We have in this paper suceeded in writing the solitary-wave solution of the SWW equa- 
tion (p]), valid for any value of the parameters a and (3, as a solitary-wave satisfying simul- 
taneously, in the slow variables, all equations of the KdV hierarchy. This is actually a result 
valid for any system presenting an exact solitary-wave solution, that is, for systems whose 
solitary-wave initial condition does not radiate. It does not matter, in this case, whether 
the system is integrable or not. 

We have then shown that, when vq is assumed to be not a solitary- wave but a two or 
more soliton solution of the KdV hierarchy, the right-hand side of Eq . (|29|) , as well as all 
other corresponding higher order equations, does not vanish anymore, which in principle 
would mean that the perturbative series does not truncate. However, in the particular 
case of a = 2/3, all those equations become homogeneous again, making then possible the 
truncation of the perturbative series through the choice of the trivial solutions 

v n = ; n = 2,3,... . (52) 

The return to the laboratory coordinates leads vq to the two-or-more soliton solution of the 
SWW equation with a = 2(3, which ultimately means that the iV-soliton solution of this 
SWW equation is nothing but a iV-soliton solution satisfying simultaneously all equations of 
the KdV hierarchy. We have in this way shown that the dynamics underlying such solutions 
can be given in terms of that governing the N-soliton solutions of the whole KdV hierarchy, 
a fact that can also be seen to be true through the AKNS scheme. 

It is important to notice the difference between the cases N = 1 and N > 2. For N = 1, 
the dynamics are equivalent for any a and f3. However, for N > 2, the above equivalence 
exists only for a = 2(3. The reason for this difference is the presence of soliton interactions 
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in the case N >2. As is well known, the SWW equation is in general non-integrable, which 
means that soliton interaction is inelastic. And it is exactly this inelastic character of the 
soliton interaction that will generate secularities other than those originated from the linear 
terms of the higher-order equations. What happens in the multiple scale method is that 
these new secularities, originated from the non-integrability of the original system, can not 
eliminated by the integrable equations of the KdV hierarchy. In this sense, the case with 
a = 2/3 is indeed peculiar, and can not be extended to the general case with any a and 
(5. In other words, being the original system non-integrable, it can not be described by the 
integrable system formed by the equations of the KdV hierarchy. 
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